Analytically integrated results for heavy fermion production in two-photon collisions 

and a high precision a s determination 
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The cross section for massive fermion production in two-photon collisions was examined at next- 
to- leading order in QCD/QED for general photon helicity. The delta function (virtual+soft) part 
of the differential cross section was analytically integrated over the final state phase space. Series 
expansions for the complete differential and total cross sections were given up to tenth order in 
the parameter /3. These were shown to be of practical use and revealed much structure. Accurate 
parametrizations of the total cross sections were given, valid up to higher energies. The above results 
were applied to top quark production in the region not too far above threshold. The cross section 
was shown to be quite sensitive to a s in the appropriate energy region. 



I. INTRODUCTION 



High energy photons may be produced by backscattering laser light off high energy e~ or e + beams. In addition, 
high degrees of polarization are possible and the photons may carry a large fraction of the electron energy. Photon- 
photon collisions also arise naturally as a background in e + e~ collisions. One major motivation for constructing a 
J> \ 77 interaction region at a high energy next linear collider (NLC) is to produce Higgs bosons on resonance via 77 
ON ■ fusion, which also allows direct determination of the f/77 coupling, which is sensitive to possible non Standard Model 
charged particles of large mass that may enter in the triangle loop. Using polarized photons allows one to control 
the backgrounds arising from 77 — > bb for an intermediate mass Higgs |Q . This background has now been studied 
' including QCD (§"P and electroweak || corrections. 
qq , In this talk we will consider in some detail the process 77 — > f f+X in the region not too far above threshold, making 
use of the complete analytical results presented in [^) , which include photon polarization. We will demonstrate the 
usefulness of 77 — » ti + X in determining a s precisely. We extend the analytical results presented in by integrating 
and obtaining analytical results for the single integral (virtual+soft) part and by series expanding the entire differential 
and integrated cross section to order /3 10 , where is the massive fermion velocity in the soft radiation limit. Such an 
expansion is shown to be of practical use, not too far above threshold, and it demonstrates many interesting features 
of the corrected cross sections. We have also provided parametrizations of the total integrated cross sections valid up 
to higher energies. 



II. GENERAL FORM AND DECOMPOSITION OF THE DIFFERENTIAL CROSS SECTION 

The process under consideration is 

7 (pi, Ax) + 7 (p 2 , A 2 ) -> f(p 3 ) + /(p 4 ) + [V(k)], (1) 

where Ai, X 2 denote helicities and the pi, k denote momenta. / (= q, I) represents a fermion with mass m and 
V = g, 7. The square brackets represent the fact that there may or may not be a gluon/photon in the final state. We 
have the following invariants, 

s=(pi+p 2 ) 2 , t = T - m 2 = (pi -p 3 ) 2 - rn 2 , u = U - to 2 = (p 2 -P3) 2 - m 2 (2) 

and 

s 2 = S 2 - to 2 = (pi + p 2 - P3) 2 - to 2 = s + t + u. (3) 

Defining 

v=l + -, w=—^- (3=^1- 4to 2 /a, x=^—^, (4) 

s s + 1 1 + p 
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we may express 

t = -s(l - v), 
s 2 = sv(l - w), 

Now introduce 



u = —svw, 

m 2 = l(i-n 



2tt- 



■7 



[Nc 



Ct = [C F ] 



ay 
2n'- 



(5) 



(6) 



where the iV c , Cp factors are present only for / = quark 
and V — gluon, respectively and e/ is the fermion's frac- 
tional charge. Here 



ay 



Oi a , V = g 
a, V = 7 



Then 



da 
dvdw 



da® daV 



dvdw dvdw 
' 1 d/(°> 



k(s) 



(i) 



27r dvdw it dvdw 



(7) 

(8) 
(9) 



The / functions are dimensionless functions of v and w, 
which allow us to parametrize our cross sections in an 
exact fashion, without dependence on ay. We use the 
normalization convention of Since, in that normal- 
ization, the /W contain an overall factor of n, we con- 
sistently present analytical results for /M /tt in order to 
cancel it. The unpolarized and polarized /W are given 
by 



f(i) 
J unp 



[/ w (+,+) + / w (+,-)], 



S=~[/ w (+^)-/ (i) (+,-)], 



in the notation /^(Ai, A2). Define 
j = 1- < AiA 2 >, 
where < A1A2 > is the average value of A1A2. Then 



(10) 



(11) 



/%') 



1+ < AiA 2 > 



1- < AxA 2 > (i) 



/ W (+,+) 
/ W (+,-) 



(12) 



= / W (+,+) + § (+,+)] 

= / (i) (+,+)-j7, 



(i) 
pol> 



so that 



/W(i) = / w (+.+) = / w (- -) 

/%•) = /«(+,-) = /«(-,+) 
/(i)( i) = / W 



(13) 



(14) 



unp 



The LO term is given by 



^4^ = * (1 -H^?(r^W (1 - 2m/ - ) (15) 



+3 



v(l — v) 



2m 2 



sv(l — v) 



5(1 -w)< 



1-/3 4 



4u 2 (l - vf 



(16) 



v(l — v) 



1-0- 



2 i 



2v(l - v) 



The last form shows explicitly the polynomial structure 
of the leading order differential cross section in terms of 
P. This is somewhat misleading, however, as we shall 
see in the next section, since the phase space in v itself 
depends on (3. 

From [^| we see that df^- 1 ' /dvdw has the form 

1 , * FJv,w) , xp/ , . . 

F h (v,w)+ ,' K ' / +F s (v)6(l-w) (17) 



tt dvdw 



Fh(v,w) 



(l-w) + 
F s (v,w)-F s (v,l) 



1 — w 

+F s (v,l)ln(l - wi)5(l - w) 
+F 5 (v)5(l-w), 



where 



w± 



4«(1 - v) 



(18) 



We have integrated analytically the terms proportional 
to the delta function part in the above equation. The 
terms Fi can be directly inferred from the expressions 
given in 0. 

It is standard B to divide the cross section (i.e. j' 1 ') 
into two parts. Firstly, there is the virtual plus soft part, 



(1) 



df (1) 

a Jy+s _ 
dvdw dvdw dvdw ' 



(19) 



where dfy /dvdw denotes the virtual contribution and 
dfcp /dvdw is obtained by integrating the bremsstrahlung 
contribution to df^ /dvdw over the region 



1 > W > Wl, so ft = 1 



sv 



(20) 



then multiplying by 5(1 — w). We follow the defini- 
tion of the soft parameter, 5, given in [|| such that the 
gluon/photon radiated becomes arbitrarily soft by mak- 
ing 5 arbitrarily small. Since this takes into account all 
virtual corrections and soft radiation, the hard radiation 
may be taken into account by integrating df^ /dvdw in 
the region w± < w < wi, so ft- Since we never reach w = 1, 
F s (v,w)/(1 — io)_|_ = F s (v,w)/(1 — w) in the hard radi- 
ation integration. We thus define df^P /dvdw as being 
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df^'/dvdw in the region wi < w < iui, so ft, where the 
(5(1 — w) terms do not contribute: 



dvdw dvdw 



(wi < W < I0l, so ft)- 



(21) 



It is not necessary to define df^ /dvdw outside that re- 
gion since it is never evaluated there. 

The sum of the (integrated over some region) hard and 
soft contributions so defined is independent of S in the 
limit 8 — > and this method of separation is referred 
to as the phase space slicing method. As one might 

expect, there is a close relation between dfy^_ s /dvdw 
and (df( 1 ^/dvdw)s as well as between df^P /dvdw and 
(df^ /dvdw)NS- Here subscript 6 denotes terms in ( p"7j ) 
multiplied by <5(1 — w) and N5 the opposite. We give 
explicitly the necessary conversion terms in Q . 

The variables v,w are suitable for performing analyti- 
cal integration of the cross section (at least for the single 
integral part). They are not, on the other hand, suitable 
for performing series expansions of the integrated cross 
section about (3 — 0. The reason is that, in these vari- 
ables, the integration limits depend on (3 so that the series 
expansion of the integrated cross section does not follow 
straightforwardly from the series expansion of the differ- 
ential cross section, and we only have complete analytical 
results for the differential cross sections. Otherwise, we 
could just expand the final integrated result. The above 
will become clear in the following sections. This approach 
also allows for cross checking; when one first expands the 
differential cross section and then integrates, the result 
should coincide with that obtained by directly expanding 
the analytically integrated cross section. We will check 
this requirement for the single integral part, for which we 
do have analytical results. 

At this point, we introduce a new set of variables, r 
and to, suitable for performing series expansions in f3. 
They are defined through 



v=-(l + r0), W =1- C/5 (t)/? 2 (1-co), (22) 



where 



We note 



1-r 2 



C ^ (T) = T-r^ 



4f (i) _ Fcpjr) ffW 
drduj 2 dvdw 



(23) 



(24) 



Because of the factor cp{r), df^/drdu; will never get a 
part proportional to 6(1 ± r). Defining 



c u = C[}(t)(1-uj)(1 + t/3), 
the invariants are given by 



(25) 



(26) 



i = -£(l_ rj 8), u = -^(1 + t/3-[3 2 c u ), 
■S2 = ^c u f3 2 , T=- S -(l-2rf3 + l3 2 ), 
U = -^[\ + 2tP + P 2 {\-2c u )], 
S 2 = S -[\ + f3 2 {2c u ~l)]. 

In terms of r and u>, the LO term has the form 

x(l + r 2 /3 2 - 2/3 2 )] (27) 
= 3(1 - w)[(2 - i)/3 + (2j + 4r 2 - 4jr 2 )/3 3 
+0((3 5 )}. (28) 

We see explicitly that the j = 0, 1 differential cross sec- 
tions, in terms of these variables, vanish by order f3 in 
the limit (3 — > 0, while the j '• = 2 differential cross section 
is order [3 3 . 



III. ANALYTIC INTEGRATION OF THE DELTA 
FUNCTION PART 

The only complete analytical results for the differential 
cross sections were presented in Analytical results 
for the virtual+soft part were presented in |^ for the 
unpolarized case and in [[| for the polarized case (where 
the virtual and soft parts are given separately, in terms 
of various functions). Still, not even the virtual+soft 
part has previously been integrated (over fermion angle) 
analytically. In this section, we present such an analytical 
integration. We were not able to integrate the non delta 
function (or hard) part analytically, in a straightforward 
fashion, and reserve that for future work. 

The integrated cross section (or /W) is obtained via 

/« = f 2 dv fdw^- = fdr f'dcj^ , (29) 
Jv! J Wl dvdw J_ x J drduj 

where 

vi = \{l-P), v 2 = \{\+P). (30) 

Let 03 be the angle between p^ and p\ in the 77 cm. 
Then #3 is given by 



cos O3 



1 — v — vw 



-V + VW) 2 + [3 2 - 1 

2t - (3c u 



^4-c u {4-(3 2 c u ) 



(31) 
(32) 



Thus, 



cos 6*3 = — - — = t, for w = oj = 1. (33) 

P 
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We see that for [3 — > 0, cos 83 varies rapidly with v, while 
it is simply equal to r. This is why the phase space in 
v becomes vanishingjy small by order (3. Similarly, from 
(^) and (^8|), or (22), we see that the w phase space is 
order /3 2 . Thus, the double integration over v and w is 
order /3 3 , in accord with (p4|). 

The integration of ( |l6| ) or (27) is rather straightfor- 
ward, yielding the LO term, 



/ (0) (i) 



2tt 
Since 



2/3(1 + /3 2 ) - 6(3 j - (1 - /3 4 + 2j) In x. (34) 



00 02/5+1 03 
l n ,= -2 E L = _ 2/3 _ 2 /L_.., (35 ) 



fc=0 



we have 

/ (0) (J) 
2tt 



(36) 



fc=2 



-1 , 1 + 2.? 



2fc-3 2fc+l 



2k+l 



so that / (0) (0,1) are order (3 and / (0) (2) is order /3 3 . 
Also, we see that only f(°'(0) is finite in the limit (3 — > 1 
and it approaches 



/C°)(0)-87T, 



for /3 -► 1 



(37) 



This is because the 1 — /3 4 term in (27) keeps the j = 
channel finite. For j = 2, the cross section vanishes for 
exactly r = ±1, as required by angular momentum con- 
servation along the 77 axis, but for j3 — > 1 the part pro- 
portional to j goes like (1 -I- t 2 )/(1 — r 2 ) as soon as we 
move away from exactly r = ±1 and is hence not inte- 
grably finite at /3 = 1. In order that the j = cross sec- 
tion be nonvanishing for r = ±1, where its maximum lies, 
the / and / must have opposite spins by angular momen- 
tum conservation, leading to m 2 / s (~ 1 — (3 2 ) suppression 
in the numerator. The fact that the LO j = cross sec- 
tion continues to be 1 — (3 2 suppressed for r 7^ ±1 follows 
from symmetry arguments [Q . This exactly compensates 
the t-channel singularity in the propagator, leading to a 
finite / (0) (0) for /3 -> 1. Of course, for r ^ ±1, the j = 
differential cross section will vanish like (1— /3 2 )/(l — r 2 ) 2 , 
making it unobservable in LO, for /3 — > 1. So, had we 
taken the limit /3 — ► 1 from the beginning, the j = 
cross section would have vanished identically. Hence the 
nonzero /(°'(0) in the (3 — > 1 limit is a remnant of using 
the fermion mass as a "regulator" . 

Near threshold, the 1 — /3 2 suppression of the j = 
channel will not be significant, hence the major con- 
straint will come from angular momentum conservation 
in the forward and backward directions which will lead 
to suppression of the j = 2 cross section there. The 
j = cross section reaches its maximum in those config- 
urations, however. Thus, we can clearly understand the 



feature of the numerical results for top quark production 
in S which show that imposing angular cuts in the di- 
rection of the beam pipe has a greater effect on the j = 
channel than on the j = 2 channel. 

We denote the single and double integral contributions 
to /« by 



(i) = 

si/di 



p V2 


dv 


f 






J W\ 


L 


dr 


f 






'0 



dw 



du> 



\dvdwj S/NS 
\dTduj J 



(38) 



(39) 



S/NS 



We performed the single integration using (j3S|), as op- 
posed to (|3^). It turned out to be quite lengthy and in- 
volved. We did not check to see whether using ([59J) sim- 
plifies the calculation. This question is probably more 
relevant to the double integration, however. Our final 
(simplified) result is 



f 1 \n(x) 
^ = am 2 /6 + a 2 Li 2 (x) + a 3 Li 2 (-x) + — 

TT (3 



a4Li 2 (a:) 



+a 5 Li 2 (-x)} + a 6 ln[(l + /3)/2][tt 2 /6 + 2Li 2 (-a?)] 
+a 7 ln(x)7r 2 /6 + a 8 {-31n 2 [(3 + /3 2 )/4] ]n(x) 
+51n[(l - /3 2 )/4] ln[(3 + /3 2 )/4] ln(a;) 
+2Li 3 [-2x/(l + /3)] - 2Li3[-2/(a;(l - (3))\ 
+2Li 2 [(l - /3) 2 /(3 + /3 2 )] ln[2/(x(l - f3))] 
-2Li 2 [(l + /3) 2 /(3 + /3 2 )] ln[2x/(l + (3)}} 
+a 9 {Li 3 [(l + /3)/2] - Li 3 [(l - (3)/2}} (40) 
+a w ln(x) + an ln 3 (x)/(3 + a±2 ln 2 (x)//3 
+ai 3 In 2 [(1 + /3)/2] ln(x) + a M p 
+ [ai 5 + a w \n 2 {x)//3 2 }(3ln{(3) 
+[ax7/3 + ais ln(x) + a w ln 2 (x)//3] ln[(l + j3)/2\. 

The coefficients Oj(j) are given in and are rational 
polynomials in j3, finite as (3 — * 0. Terms of the form 
F(f3) - F{-p) will vanish for /3 = 0. Noting that 



Li 2 (l) 



Lia(-l) 



12' 



(41) 



we see that only the terms proportional to ai, . . . , 05 
may contribute at threshold, which is the case for j = 0. 
Indeed, one finds the correct threshold correction from 
those terms alone. The relevant series expansions will be 
given in the next section. As we shall see in Section V, 
the double integral series expansion starts at order /3 3 . 

Two independent determinations of (|4C| ) were made us- 
ing Mathematica @ and REDUCE jnf. That software 
could not evaluate certain integrals which can be found 
in jl2| . It was verified that the analytically integrated re- 
sult agreed numerically with the numerically integrated 
result. In the next section, we will show how one can use 
the series expansion as a very solid check as well. 

Perhaps the most convincing check of (^) and the 
analytical result for df^/dvdw (or dajsiho/dvdw) ob- 
tained in p| is the excellent numerical agreement with 
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tabulated results for /W existing in the literature. The 
only existing analytical results, aside from those in [Q, 
are the expressions for df^l 3 /dvdw (i.e. /dvdw) 5 ), 

dfg /dvdw given in Q for the unpolarized case (us- 
ing dimensional regularization), with which we agree 
exactly, and similar expressions for the polarized case 
in H (obtained using a gluon energy cut and a small 
gluon mass as infrared regulator). The latter are not 
quite in a form suitable for direct analytical comparison. 
There have been no other analytical results presented 
for (df^ I 'dvdw)N& in the polarized or unpolarized cases. 
Hence we must perform the above mentioned numerical 
checks. 
Define 



with @ for /W (+,-), we conclude that our analytical 



results for df^/dvdw of [|2j have been verified. In light 
of the above, Table || is seen to be the most complete and 
precise such table at present. 



2m 



1 



(3 = y/l- 1/z 2 



(42) 



(i) 



In Table [j] we give numerically computed values for / U np, 
/poi' /^(+> +)> ~) as wen as tne specific contri- 

butions from all the and fj^ to the corresponding 
/W, for various values of 1.2 < z < 20. The result at 
z = 1 is given exactly by the series expansions presented 
in the next section. We also indicate the number of sig- 
nificant figures, n.s., following the decimal point, in fjj, 
(and /«). 

The next issue is, of course, how well these values com- 
pare with other tabulated values for / W . Two other such 
tables exist at present. The original one of [Q gave / unp 
for z = 2, 3, 4, 5, 10; the value at z = 1 being numerically 
equal to the known threshold result, as given in the next 
section. Their numerical values were obtained using the 
/v+'s" 1 ' gi ven m [||: added numerically to /F np , deter- 
mined there using the same methodology as |g] , which is 
equivalent to our method. We find numerical agreement 
with to within the precision of those values, which is 
roughly at the order of one part in 10, 000 or better. This 
can only be achieved with correct analytical results. Our 
calculation of is identical in method (same integrals 

and structure) to that of / unp (at the differential and in- 
tegrated level), the only difference arising from different 
traces due to the contraction with a polarized photonic 
tensor rather than an unpolarized one. As two indepen- 
dent determinations of these traces were performed, there 
is little room for any error in Fortunately, we may 

directly check this assertion since the values of /^ (+, +) 
and / (1) (+, -) for z = 2, 3, 4, 5, 10, 20, 50 were tabulated 
in ||. There, Monte Carlo methods were used, leading to 
accuracy at the level of better than 1% in regions where 
the /W are sizable, but apparently not better than ±0.2 
or so in absolute error. This absolute error is noticeable 
only for /^ (+, +) and only for z — 2, 3, where /W (+, +) 
is small. To within the above accuracy, we are in good 
agreement with j|. Since f (1 \+,-) = /inp- fpj and 
since we have precision agreement with for / un p and 



IV. SERIES EXPANSION OF THE DELTA 
FUNCTION PART 

Besides providing a useful check of the analytical inte- 
gration of the previous section, there are many reasons 
why it is useful and instructive to series expand the dif- 
ferential and integrated cross sections about (3 = 0. In 
the absence of complete analytically integrated results, 
only a series expansion about f3 = can be used to make 
(very) high precision predictions in the (3 ~ region. 
One also sees the structure of the cross section in a way 
that cannot be inferred from the non-expanded analytical 
results, which are somewhat complicated. From a practi- 
cal viewpoint, having "simple" series expansions for the 
differential cross sections allows one to do complete nu- 
merical studies in the region not too far above threshold 
rather easily. This is because the resulting expansions 
only involve simple polynomials and simple logarithms. 
We will address the issue of the region of validity of the 
expansions as well. 

The other issue is that of resummation. There are large 
correction terms at threshold which can be resummed. 
Having a series expansion of high enough order to be of 
practical use allows one to explicitly perform resumma- 
tions up to some order in f3 while leaving the higher order 
terms the same. The net result would be an equally sim- 
ple series, improved via resummation so as to allow one 
to go closer to threshold. This is beyond the scope of 
this paper as are other very near threshold effects. Suf- 
fice it to say that having the threshold series expansion 
will facilitate these studies for those interested. 

Throughout, we will expand up to order /3 10 (includ- 
ing j3 11 \n(3 terms). The expansion which exists in the 
literature (see 0) is only for / un p and only goes to or- 
der f3. Going to order (3 W may seem excessive at first, 
but we found it to be a good stopping point for several 
reasons. Considerable structure arises beyond order /? 
which allows us to see the general, all-orders in f3, struc- 
ture of the various series. Also, one gains little in terms 
of precision by going to even higher orders in (3, with- 
out including several more terms. Then, the series would 
start to become lengthy and cumbersome, reducing the 
advantage over the analytical result in terms of ease of 
use. For certain series, going much beyond (3 W would 
take a very large amount of computer memory and run- 
time, not justifying the extra effort, as going to order /3 10 
was a considerable task in itself. Finally, by going to such 
a high order, we may stringently check the analytically 
integrated single integral result of the previous section as 
will be described below. 
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We find that d/W /drdu may be expanded in the gen- 
eral form 



drdoj 



(43) 



i=0 i=0 



Therefore may be expanded as 



f(i) 



(44) 



»=o i=o 



where the dij are given by 



dij = dr duo Cij(r,u)). 
J -I Jo 



(45) 



With the variables v and w, the integration limits depend 
on 0, hence the above arguments do not hold. So, one 
sees clearly the necessity of the change of variables. 

We convert from (d/W / 'dvdw)$ to (dp 1 ) / 'drdu>)s using 
(E4j), which modifies the overall factor via 6(1 — w) — > 
0^cp(r)6(l - w)/2 = 06(1 - u)/2. The compact results 
for the series expansions of (d/W / 'drduj)s for j = and 
j = 2 can be found in || , up to order W . 

We notice that the cross section is isotropic up to or- 
der 0; the angular (t) dependence enters only at order 
2 . The LO term, on the other hand, was isotropic up 
to order 2 . We also see that the step function threshold 
behaviour arises entirely from the j — channel, at the 
level of the differential cross section, since the j = 2 chan- 
nel starts at order 2 . From the 1 — r 2 overall factor, and 



using (|33|) , we see that the delta function contribution to 
the j = 2 cross section vanishes at cos #3 = ±1 as did 
the LO cross section (|27]). This vanishing is not obvious 
from the exact analytical expressions, but simply reflects 
angular momentum conservation along the 77 axis when 
Lu = 1 (2 — > 2 kinematics). The j = channel, on the 
other hand, becomes infinite (but integrably finite) for 
cos 83 = ±1 due to the ln(l — r 2 ) terms. 

The expansions have rather simple structure in that 
the Cij(r) are simply polynomials in r. This amounts 
to considerable simplification and reduction in computa- 
tional time relative to the exact expressions, especially 
after the non delta function part is added, where the 
simplification is even greater as we shall see in the next 
section. 

Two independent calculations of series expansions were 
performed using Mathematica and REDUCE. The ex- 
pansions were also checked numerically by subtracting 
them from the exact expressions. The difference was 
checked to be of order 11 . This is most straightfor- 
wardly done by taking rather small 0. 

Assuming we are working at where the series are 
sufficiently accurate, one could easily analytically inte- 
grate them over a region of r (cos #3) relevant to some 
experiment, if desired, and implement angular cuts ana- 
lytically. After a suitable change of variables, the same 



could be done for the hard radiation part, either ana- 
lytically or numerically. Cuts on additional observables 
may be made by subtracting off the unwanted configura- 
tions using the squared amplitudes given in Q and Monte 
Carlo integration, for instance. Here, we simply present 
the total integrated results. 
For the j = channel, we find 

k ( t x) (+,+) = 

2 j2vr 2 - (20 - tt 2 )0 + 10/3vr 2 /3 2 + 3 /3[- 340/3 + n 2 

+641n(2/?)] + 8/f57r 2 /3 4 + 4/1 5/3 5 [-6343/45 - tt 2 
-321n(2) + 256/3 ln(2/3)] - 104/105tt 2 /3 6 
+4/105/3 7 [-39163/315 - vr 2 + 208/3 ln(2/3)] 
-88/3157T 2 /? 8 + 4/9/3 9 {l/7[-7r 2 /5 - 64/3 ln(2)] 
+ 1/25[128 \n(20) - 43903/315]} - 488/3465tt 2 /3 10 



-103232/51975/3 n ln(2/3) 



(46) 



(47) 



and for j = 2, 
(+,") = 

TT 

16/3|tt 2 /3 2 - 80 3 + tt 2 a + 32/3 5 [-289/720 + ln(2)/5 

+ ln(2/3)/3] + tt 2 /70 6 + 6/50 7 [-3947/945+ 16/7 ln(2) 
+32/9 \n(20)\ + 29/105tt 2 /3 8 + 4/15/3 9 [-823/45 + 8 ln(2) 

+16 ln(20)] + 289/11557T 2 / 3 10 + 256/63/? 11 ln(20) 



The results are indeed quite simple. We may again obtain 
^/v+s by adding the corresponding conversion term ||. 

The strongest check comes from the fact that the ex- 
pansions ((46|), ( [i"7| ) which come from integrating (|4^), 
( pli] ) agree exactly with the expression obtained by ex- 
panding the analytically integrated result ( f4(i| ) directly. 
In this way, we simultaneously check all the above men- 
tioned expressions, including our analytical integration 
©. The expansions (f46|), ( [47| ) were also checked nu- 
merically by subtracting them from the exact expression 
( pp| ) and verifying that the difference was order 11 . 



V. SERIES EXPANSION OF THE NON DELTA 
FUNCTION PART 

Perhaps the most remarkable result of the series expan- 
sion is the simplification of the non delta function part, 
whose original form is the most lengthy part of the exact 
result, involving complicated logarithms, etc. . . Although 
the intermediate expressions were very lengthy and con- 
siderable computational time was required, a large degree 
of cancellation resulted in the following simple series. We 
convert from (df^- 1 ' /dvdw)isrs to (df^ 1 * /drdw)jvs by mul- 
tiplying by 3 C/3(t)/2. Again, we present here results 
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for the total cross sections, the results for the differential 
cross sections can be found in |J. Again, two indepen- 
dent determinations were performed using Mathematica 
and REDUCE. The expressions were also checked numer- 
ically analogously to the delta function part. The inte- 
gration of differential cross sections over r, u> is straight- 
forward and we obtain 



34624 
2205 



42368 



3675 



-1024 
~45~ 



2816 
525 



P 7 



134656 
19845 



(48) 



f3 9 . (49) 



We notice the absence of any logarithms, including 
powers of In (3. The structure is fairly predictable as well. 
We see that the series begin at order /3 3 and (3 5 for j = 
and j — 2 respectively, so that their effect will be negligi- 
ble very near to threshold. On the other hand, the large 
coefficients imply that they soon become noticeable for 
small (3. 

These are remarkably simple results, which suggest 
that the exact integrated result for fjj^ is not too compli- 
cated. We notice the vanishing of the coefficients of the 
even powers of (3. This follows from the antisymmetry 
in r of the corresponding terms in the differential cross 
section. For discussion about the conversion terms see 



VI. TOTAL SERIES RESULTS AND 
NUMERICAL PARAMETRIZATIONS 

We are now in a position to study the total cross sec- 
tion, by combining the results of the previous sections. 
Adding (^) and ( [48"1 ) gives the series for the j = total 
cross section 



-/ (1) (+,+) = 

7T 



2 |2tt 2 - (20 - tt 2 )/? + 10/3tt 2 /3 2 + ^/3[-404/3 + tt 2 

+641n(2/3)] + 8/15tt 2 /3 4 + 4/15/3 5 [-6511/45 - tt 2 
-32 ln(2) + 256/3 ln(2/?)] - 104/105tt 2 ^ 6 
+4/105/? 7 [25757/315 - tt 2 + 208/3 ln(2/3)] - 88/315tt 2 ^ 8 
+4/9/3 9 {1/7[-tt 2 /5 - 64/3 ln(2)] + l/25[1281n(2/3) 
+407639/2205]} - 488/3465tt 2 /3 10 



-103232/51975/3 n ln(2/3) 



(50) 



and adding (^7|), (^ 
cross section 



gives the series for the j = 2 total 



-/ (1) (+,-) 



16/3|tt 2 /3 2 



8/T + tt 2 /? 4 + 32/3 5 [-77/144 + ln(2)/5 



+ ln(2/3)/3] +tt 2 /7/3 6 + 6/5/3 7 [-677/135+ 16/71n(2) 
+32/91n(2/3)] + 29/105tt 2 /? 8 + 4/15/3 9 [-16949/735 
+8 ln(2) + 16 ln(2/?)] + 289/11557T 2 /? 10 

+256/63/3 11 ln(2/3) 



(51) 



Such simple expressions indeed make numerical studies 
not too far above threshold rather straightforward. We 
can get an idea of how well these series work for typical (3 
by comparing with numerically calculated values of f^ 1 '. 

In Table || we present the fractional error on the se- 
ries for / (1) (+, +), / (1) (+, -), /un P relative to the result 
obtained using numerical integration, for various values 
of z in the region 1.05 < z < 1.4. For z < 1.05, the 
series expansions are more accurate than the numerical 
results. At z = 1.05, the errors are at the 10 -7 — 10~ 6 
level. For z = 1.2 they are at the 10 -4 — 10~ 3 level and 
for z — 1.4 they are at the 10 -3 — 10~ 2 level. The errors 
on f^ 1 '{+,+) are at the lower end, while the errors on 
/^(+, — ) are at the higher end and those for lie 
in between. This is good because, as we shall see in the 
next section, in determining ot s via top quark production 
at a 77 collider, it is the j = and unpolarized channels 
which are of interest, the j — channel being the most 
interesting one. With precision of better than one per- 
cent for z < 1.4, we have sufficient accuracy to use the 
series expansions (differential in particular) to perform 
easy numerical studies relevant to top quark production 
at a 77 collider of y/s < 500 GeV. As we shall see, for the 
a s determination, going to much higher energies is not 
useful since the determination is best done near z = 1.2 
(V5~ 420 GeV). 

It is also useful to be able to parametrize to good 
accuracy for larger (3, relevant for bottom and charm 
quark production at intermediate energies or top quark 
production at very high energies. This was done by fit- 
ting numerically computed values of /W. We divide 
the parametrizations into 3 regions: a low energy re- 
gion (1 < z < 1.5 or < f3 < 0.7454), an intermediate 
energy region (1.5 < z < 5) and a high energy region 
(5 < z < 20). We will denote the corresponding as 
/W' ie , f^- le and /(!)^ e , respectively. 

The various forms for the parametrizations are 



/ (1) '' e (+,+) = 27T 



2^ 2 



(20-tt 2 )/3 + 



10tt' 



ft 2 



-^Vln/? 
3 



7 

£ 

i=3 



5>(*-1.5)*, 
4 



(52) 



i=0 
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and 

ifi 10 

/ (1) < Ze (+, -) = if [- 2 /3 2 - 8/3 3 + tt 2 /? 4 ] + £ c#, 

4 

/ (1) ' <e (+,-) = 1.5)', 

i=0 

-) = £ Ci ( z -5)\ (53) 

i=0 

The Ci are given in J^]. In the low energy region, where 
high accuracy is required, the parametrizations are accu- 
rate to < 0.01%, with the errors being the largest near 
the higher end of the region. The leading terms, given 
analytically guarantee the correct threshold behaviour 
as they are just those in the exact series expansion. As 
mentioned earlier in connection with the series expan- 
sions, one can explicitly perform resummations on those 
terms. Thus one could modify the above parametriza- 
tions to include resummation effects without changing 
the higher order coefficients. 

In the intermediate energy region, f^' le {+, — ) is accu- 
rate to < 0.1%. / (1) ' 4e (+, +) is accurate to < 1%, except 
very near / ( - 1 - ),le (+, +) = 0, which occurs for z ~ 2.15, 
3.15. There, the absolute errors remain small, but of 
course the fractional error is larger. In the high en- 
ergy region, f^ 1 ' ,he (+,— ) is accurate to < 0.05%, while 
/ (1) ' ,le (+,+) is accurate to < 0.5%. The above errors 
are rather conservative and one can not distinguish the 
parametrizations from the exact results for practical pur- 
poses. 

VII. PRECISION as DETERMINATION FROM 
TOP-QUARK PRODUCTION 

A high energy 77 collider can be used as a "factory" 
for many interesting particles: Higgs bosons, W bosons, 
top quarks etc... The beam polarization we be useful in 
producing Higgs bosons and reducing QQ backgrounds. 
More specifically, the j = channel will be of interest. 
This channel also turns out to be the channel of interest 
when trying to determine a s via top quark production, 
making it complementary to the Higgs studies. The rea- 
son is that the cross section, and QCD corrections, are 
enhanced in this channel, thereby improving the statis- 
tics and the determination of ct s , to which the cross sec- 
tion will be quite sensitive. The process 77 — > tt + X 
is more powerful than e + e~ — * QQ + X in determining 
a s because the QCD corrections are quite small in the 
latter, thus requiring an unreasonably large number of 
events for high precision; the corrections are suppressed 
by a s /ir ~ 4%, relative to the Born term. In 77 — > ti+X, 
we can "pick" our QCD correction by choosing the ap- 
propriate beam energy. Of course, as one gets too close 
to threshold, the perturbation series cannot be trusted, 



for reasons we will discuss below. Hence there are limi- 
tations. 

To best illustrate the above idea, in Fig. [I] we have 
plotted the 77 — > tt + X cross section at LO and NLO, 
in the region 1 < z < 1.4, for the various helicity states. 
We took N f = 5 , m t = 174 GeV and used A = 230 MeV 
in the two- loop expression for a s , evaluated at [i 2 = s. 
One could also use Nf = 6, but since we are not far above 
threshold it is simpler to use Nf — 5 for evolution from 
/i 2 = Af| to /1 2 — s. We notice that the j = cross 
section is the largest, as are its QCD corrections, in this 
region. The region z ~ 1.2 is nice in that the j = cross 
section is near its maximum and the QCD corrections 
are sizable (~ 20% of the total cross section), yet not so 
large that the perturbative expansion is unreliable. As 
one gets closer to threshold, other higher order effects, 
nonperturbative effects and top width effects may also 
become important. For these, and other reasons to be 
considered below, we will suggest z — 1.2 as being the 
optimal region for extracting a s and we will give a rough 
estimate of how precisely a s may be determined there. 
As well, we suggest the j = channel as being the most 
powerful. 

Firstly, we note that z = 1.2 corresponds to ^/s 77 ~ 
420 GeV, for top quark production. This energy should 
be accessible at a y/s e + e - > 500 GeV NLC. A typical 77 
luminosity assumed is 20 fb _1 . Since a ~ 1.4 pb, this 
corresponds to roughly 28,000 tt events. Since the QCD 
correction is ~ 20% of the total cross section, this trans- 
lates to Aa s /a s ~ 3%, statistically. With a luminosity 
increase and, possibly, extended running, one could en- 
vision going to the percent level or better. 

The above analysis was purely based on statistics and 
one-loop QCD corrections. Therefore, we will briefly dis- 
cuss various theoretical systematic uncertainties. Clearly 
one needs a two-loop analysis when dealing with one-loop 
corrections of order 20%, in order to determine a s at the 
level of a few percent. Threshold resummation can also 
be performed. One should also take into account the 
one-loop electroweak corrections ||. The QED ones are 
identical in form to the QCD ones, with the appropri- 
ate change in normalization, given by (Q) . There will be 
a minor dependence on m t , which will be lessened with 
future Fermilab runs. The uncertainty on m t translates 
to an uncertainty on z. Since the j = cross section is 
near its peak for z ~ 1.2, minor variations in z will not 
appreciably affect the results. 

Of some concern are resolved photon contributions, 
where a gluon or quark within the photon can participate 
directly in the interaction. Suppression of these contri- 
butions is a major reason for working close to thresh- 
old. Since the parton distributions within the photon 
drop steeply with increasing momentum fraction, a;, and 
since x must be large near threshold, such contributions 
are quite suppressed. Confirmation of this assertion may 
be inferred from the resolved contributions to b quark 
production near threshold presented in Jl3| from which 
we conclude that only very poor knowledge (if any) of 



8 



the photon structure will be required, as such contribu- 
tions will be a fraction of a percent of the cross section. 
One can further reduce those contributions by identify- 
ing outgoing jets collinear with one of the photon beams, 
which are a signature of resolved photon events. One 
can also require that the energy deposited in the detec- 
tors be equal to the total beam energy in order to account 
for missed jets of the type mentioned above. 

From the experimental side, we are assuming only 
that tt events can be clearly identified. With experience 
gained from Fermilab, this seems reasonable, especially 
considering the cleaner initial and final states in the 77 
case. Another experimental issue is that of normaliza- 
tion. We suggest the measurement of a ratio 



W1 _ (7(77 ^QQ + X) 



(7(77 -> PP + X) 



P = W,l. (54) 



The ratio of tt to W + W~ events is statistically quite 
powerful as over one million W + W~ events are expected 
at such a "W factory" (b|. This highlights the comple- 
mentary nature of top quark and production at a 
77 collider. As well, electroweak corrections to W + W~ 
production have been studied [M . For the same reasons 
as for tt production, the resolved photon contributions 
will be suppressed. If a bb pair is produced in conjunc- 
tion with the W + W~ , this will constitute a background 
to tt production. 

It is worth discussing the many advantages of deter- 
mining a s via 77 — > tt+X relative to some of the options 
currently being used. The calculation is perturbative and 
avoids nonperturbative contributions arising in a s deter- 
minations from mass splittings and tau decays. Other 
determinations, based on evolution of hadronic struc- 
ture functions, rely on the parton model and assumed 
knowledge of hadronic structure. No such assumptions 
are made here. Unlike the 3- to 2-jet ratio from e + e~ 
annihilation, we avoid having to define the jet isolation 
criteria by measuring the total tt cross section. Since we 
are at a large energy scale, not only does perturbation 
theory work well, but we automatically determine a s at 
(or above) the tt threshold, without having to perform 
evolution or cross flavor thresholds. From a theoreti- 
cal viewpoint, the most comparably clean determination 
comes from the ratio of hadrons to lepton pairs produced 
in e + e~ annihilation at the Z pole. As mentioned ear- 
lier, the small QCD correction proves an insurmountable 
limiting factor in that case. 

At this stage, our enthusiasm is dampened somewhat 
however by the need for a two-loop calculation. This 
need is highlighted by the fact that there is an arbitrari- 
ness in the choice of renormalization scale, fi, which can 
only be compensated by the inclusion of two-loop correc- 
tions. The variation of a s with ln/x is order a 2 s though, 
so for a reasonable choice of \i (i.e. y/s, mt, • ■ ■ ) the 
two-loop scale dependent contribution should not be too 
large and should not change the value of a s radically. 
Nonetheless, as pointed out earlier, a two-loop calcula- 



tion will eventually be required. In light of that fact, we 
see the importance of having simple analytical results for 
the one-loop corrections as they will be incorporated in 
the two-loop result. Also, we do not suggest that one 
consider this determination of a s in isolation. Rather, it 
should be combined with all other precision determina- 
tions, including the low energy ones, in order to minimize 
the error and provide an excellent test of QCD. 
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FIG. 1. The 77 -» tt + X cross section at LO and NLO, 
versus z, for the various helicity states. 







TABLE I. The various for values of 1.2 < z < 20, and the corresponding single and double integral contributions. Here 
n.s. is the number of significant figures after the decimal point in and f£\ 



z 


n.s. 


/£> ('+>+) 




At) 

J unp 






/i 1)d (+°-) 


/«(+!-) 


1.2 


4 
4 


70.578894 
103.9951788 


-5.47998017 
-6.85659845 


65.0989 
97.1386 




33.4162848 
37.1626092 


-1.37661828 
-4.10336189 


32.0397 
33.0592 


2 


3 
3 


68.5516 
-7.8931 


-24.064 
14.728 


44.488 
6.835 




-76.4447 
144.9963 


38.792 
-62.856 


-37.653 
82.140 


3 


3 
3 


92.2075 
-99.5479 


-29.594 
96.259 


62.614 
-3.289 




-191.7554 
283.9629 


125.8526 
-155.447 


-65.903 
128.516 


4 


3 
3 


132.0495 
-153.7108 


-33.0381 
182.4102 


99.011 

28.699 




-285.7603 
417.8098 


215.4483 
-248.4864 


-70.312 
169.323 


o 


Q 
O 

3 


1 7fi 701 A 

-190.8526 


-ou.ouz 

263.100 


72.247 




Q«7 KKAf\ 
-OU I .OOIU 

544.2554 


9QQ Q09 

-336.704 


-u t .\jo£ 
207.551 


10 


3 
3 


395.3262 
-292.8618 


-55.3625 
584.182 


339.964 
291.320 




-688.1880 
1083.5142 


639.5445 
-694.907 


-48.6435 
388.607 


20 


1 
1 


749.8886 
-390.5080 


-79.437 
1007.530 


670.45 
617.02 




-1140.3966 
1890.2852 


1086.967 
-1166.404 


-53.43 
723.88 


TABLE II. The fractional errors 
1.05 < z < 1.4. 


on the various 


computed using 


the 


series expansions up to order f3 10 , 


for values of 


z 




P 


f.err(+,+) 


f. err(+, 


") 




f . 6IT un p 




1.05 
1.1 
1.2 
1.3 
1.4 




.3049 
.4166 
.5528 
.6390 
.6999 


2.1 x 10~ 7 
-6.8 x 10~ 6 
-2.0 x 10~ 4 
-1.4 x 10~ 3 
-5.6 x 10~ 3 


2.5 x 10 
2.3 x 10 
3.3 x 10 
1.3 x 10 
2.9 x 10 


-6 
-4 
-3 
-2 
-2 




4.3 x 10~ 7 
3.1 x 10~ 5 
6.9 x 10~ 4 

3.4 x 10~ 3 
8.9 x 10~ 3 


2.1 x 10~ 6 
6.6 x 10~ 5 
1.5 x 10~ 3 
7.3 x 10~ 3 
2.0 x 10~ 2 
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